UND-HEP-OO-BIG 03 
hep-ph/0005278 



Analytical 



Heavy Quark Expansion in the 



't Hooft Model 



Matthias 



Burkardt a and Nikolai Uraltsev & ' c 



a Dept. of Physics, New Mexico State Univ., Las Cruces, NM 88003-0001, U.S.A. 
h Physics Dept., Univ. of Notre Dame du Lac, Notre Dame, IN 46556, U.S.A. 
c Petersburg Nuclear Physics Inst., Gatchina, St. Petersburg, 188350, Russia 



We present a number of exact relations for the heavy quark limit and develop an 
analytical 1/tuq expansion for heavy mesons in the 't Hooft model. Among the new 
results are relation = A — m^ p + (3 2 , 1/raq corrections to the decay constants, to 
the kinetic expectation values and l/mg nonperturbative corrections to the B — ► D 
amplitude at zero recoil. The properties of the IW functions are addressed and the 
small velocity sum rules are verified. 



Abstract 



1 Introduction 



Heavy quark symmetry and the heavy quark expansion have played an important 
role in understanding weak decays of charm and beauty hadrons allowing extraction 
of the fundamental parameters of the Standard Model. Beauty and, in particular, 
charm quarks are not infinitely heavy even in the crude approximation. In practice, 
1/mg corrections to the strict tuq — > oo limit often constitute the main limitation. 
Even in a few cases where the first terms are known, a question remains about the 
convergence of the employed l/mq expansion. (For a review and further references, 
see Ref. |TJ.) 

In such a situation it is advantageous to have a model laboratory where both l/mq 
corrections, and the whole finite-mg amplitudes can be evaluated exactly. On the 
one hand, this allows to trace in detail how the methods used for actual QCD work in 
a simplified setting. On the other hand, it has been empirically observed that certain 
quantities in the heavy flavor hadrons suffer from numerically large power corrections, 
whilst others seem robust against finite mass effects. Studying this in the toy models 
can help gaining some insights which can be applied in various studies of actual QCD, 
including those based on the QCD sum rules technique or lattice computations. 

One such solvable model, which has been applied to a variety of strong interaction 
phenomena in the past is the 't Hooft model, QCD in 1+1 dimensions at large 
number of colors iV c — > oo [0]. The 't Hooft model has two important features 
that resemble real QCD 3+1 phenomenologically: confinement and spontaneous chiral 
symmetry breaking. Since the underlying microscopic mechanisms for these effects 
in this model are quite different from the respective mechanisms in QCD3 + i, it is 
probable that QCDi + i is of little help to understand the origin of these features 
in real QCD, but it may nevertheless be very valuable to better understand the 
consequences for other observables. 

In this paper we will focus on the tuq — > oo limit of the 't Hooft model. Even 
though we were not able to solve this limit analytically, we succeeded to derive a 
number of exact relations involving terms that appear in the 1/rriQ expansion. In 
Section 2 we will discuss some general features of this model. In Section 3 the static 
limit rriQ — > oo will be explored and we will derive a number of relations involving 
't Hooft wavefunctions in this limit. These relations will be applied in Section 4 to 
derive some exact results for power corrections in this model. Section 5 gives the 
conclusions. 

2 The 't Hooft Model 

QCDi + i is based on the Lagrangian 

£ = £g(*7 M A*-™ 9 )9 " ^TrG^GV, (1) 

q ^9s 
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defined in one space and one time dimension, i.e. //, v = 0, 1 and all fields depend 
only on two space-time coordinates. In the light-cone gauge, A + = 0, as in any other 
axial gauge, the non-Abelian term in the field strength tensor = i [_D M , D u ] = 
d^Ay — d u A^ + i [Afj,, A u ] vanishes, which drastically simplifies the dynamics of the 
model. This allows one to eliminate the only non-vanishing component of the gauge 
field A~ , by means of the Poisson equation 



d 2 _A- 



(2) 



where we have introduced light-cone coordinates x ± = 4= (x°±x 1 ) . In two dimensions 
there are no dynamical gluons and, after solving the Poisson equation (0), the only 
remnant of the gluon field is a Coulomb like instantaneous interaction among the 
quarks. Of course, in one space dimension, the 'Coulomb' interaction is linearly 
confining. It is this feature which makes QCDi + i so attractive if one is interested in 
studying models which exhibit infrared slavery. 

An additional simplification occurs in the large N c limit 



Nr. 



gjNc 

2tt 



fixed, 



(3) 



where sea quarks (more precisely, quark loops) are suppressed and only planar di- 
agrams survive. Therefore, in a light-cone Fock space expansion, the valence qq 
approximation for mesons becomes exact, and the two body equation for the quark 
distribution amplitude <f n ( x ) of mesons (the 't Hooft equation H) reads 



m\ — (3 2 m| — (3 2 

x 1 — x 

m\ — [3 2 m\ — f3 2 



(p n (x) - (3 dy 
Jo 



x 



1 — X 



(y - x) 2 

(f n (x) + V if n {x) . 



(4) 



Here x denotes the light-cone momentum fraction carried by the quark and M n is the 
invariant mass of the meson. 

Eq. (!]) has a very physical interpretation: the light-cone energy of the qq pair 
consists of a sum of the light-cone kinetic energies of the quark and antiquark plus 
an interaction term V. The integral operator V can be interpreted as the mo- 
mentum space representation of a linear potential. The singularity of the QCD- 
Coulomb interaction in Eq. (f|) is regularized using the principal value prescription, 

with/ 1 dy^5 I = -(i + 3^). 

For practical purposes, it is convenient to rewrite Eq. (U) into the form where the 
singularity of the interaction term is less severe: 



* = — + 

V x 1—x 



ip n {x) - f3 2 



1 dy ¥n{y)-¥n{x) 

o (y-x) 2 



(5) 
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Note that Eq. (|5|) can also be obtained as a variational equation for the functional 
(Hamiltonian) defined by the quadratic form 

The interaction term in Eq. (H) is non-negative and the kinetic term is minimized for 



mi 



rl _ n — ^, i.e. when the ratio of momentum fractions carried by the quark and 
antiquark equals their ratio of bare masses. Therefore the variational formulation 
Eq. (H) of the 't Hooft equation yields a lower bound 

M n > m l + m 2 (7) 

(the inequality is saturated only for the massless pion where m\ = m<i = and at 
n = 0, with (p^(x) = 1). The mass of a bound state exceeds the sum of the bare 
masses of the constituents, the fact expected semiclassically. 

A certain parity relation proves to be useful, it was first given in Ref. |J and 
reportedly ascends to 't Hooft. We define operator K as 



Kip(x) = [ dy 
Jo 



y-x 

Then the following commutation relation holds: 

y 1—x Jo L—y 



x Jo y 1—x Jo l — y 



Since for eigenstates the expectation value of the commutator of any operator A with 
Hamiltonian vanishes, 

(n\[H,A]\n) = 0, (10) 

one has 

= (n\[H,K]\n) = m \(j^ Vn { y )\ - m \ (£ ^- Mv)) ■ (11) 



Therefore, 



m 2 I -p-Vniy) = -P n m x [ —<p n (y), (12) 
Jo l—y Jo y 

and P n is identified with the parity of the eigenstate n. This identification is confirmed 
by comparing with explicit expressions for the matrix elements of the pseudoscalar 
and scalar densities between the state n and the vacuum 101 



(0\mi\n) = I dx(--—) 



\q2ij5qi\n) = J o dx + ) <p n (x). (13) 
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Using Eq. (|12"D , it is easy to see that (0|g2<7i|^) 7^ only for states with P n = 1 and 
(0 1 g2^75<?i I n) 7^ only for states with P n = — 1. 

Integrating the 't Hooft equation over x and using the parity relation (12) we get 

Ml n 



1 dx 



X 

1 da; 



mAmi- P n m 2 ) Jo 



Ml 



dx ip n (x) 



dx (f n {x) 



lo l — x r "' y ~ / m 2 (m 2 — P n mi) Jo 
Similarly, integrating the 't Hooft equation multiplied by x yields 

m 2 



(14) 



dx x(p n (x) 



dx (1 — x)ip n (x) 



1 1 
m\ — vri2 



M 2 

n 

1 1 

m 2 —m\ 
Ml 



+ 



+ 



m 2 -P n m 1 

Uli 

m 1 -P n m 2 



[ dxLp n (x)--^— [ dx \n- — -ip n {x) 
Jo Ml Jo x 



dxtp n (x) + 



Ml 







1 l — x 

dx In (p n (x)- 

x 



Another useful relation for the 't Hooft equation which holds for arbitrary quark 
masses was derived by Burkardt M and is often referred to as the virial theorem. Let 
us denote 

D = 4- b = ^w=d ■ < 15 » 

A direct computation yields for the commutator of D with the interaction operator 
V in the 't Hooft equation the following result: 

[D,V] = -V (16) 

and, therefore, 

[D,H] = H ■ \ (17) 

Likewise 

[D,n] = -n + 



ml -ft 2 
(1-x) 2 

mj-p 2 

x 2 



Using again the fact that the expectation value of the commutator of any operator 
with Hamiltonian vanishes in an energy eigenstate, these two commutation relations 
yield the relation which holds for any eigenfunction:0 



Ml 



K-/? 2 )<777> = K-/? 2 )( 



(1-x): 



(19) 



The 't Hooft equations describing mesons with one heavy quark Q acquire new 
symmetry which entails one to a number of additional relations. In the rest of the 
paper we study these relations in the static limit (jtiq — > 00) and develop the 1/mg 
expansion around this limit. 



1 As it stands, Eq. ( |l9| ) is valid only for m 2 > (3 2 . For m 2 < (3 2 it can be replaced by a subtracted 
version j|. 
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3 Static version of the 't Hooft equation 

In the heavy quark limit, i.e. when the mass of the quark mi — ► oo, the meson wave 
functions become sharply peaked near x — > 1, since most of the momentum is then 
carried by this quark. Therefore, in order to study the mi = itlq — > oo limit, it is 
convenient to introduce the nonrelativistic variables M n = mQ + e n , t = (l-i)mg and 
ty n (t) = ^=(p n (l — -^^j. In these variables the 't Hooft equation <Q) takes the form 



2m a 




a 2 



as 



(t-s) 



■ (20) 



We assume that in the notations of the previous section mi = mQ is large, and 
m sp = ni2 will be denoted simply as m below. The limit mQ 3> (3 is obtained by 
expanding the second term in the r.h.s. of Eq. (p0[) in t/mq and extending it to the 
interval [0, oo) Jl 



e n + 



2m a 



2t 



t t 2 t 3 
2 + 2m7 ! 



2m^ 



^ 2 



00 , tf„(s) 
ds 



(f - S ) 



Performing 1/mQ expansion, it is convenient to study the eigenvalues 



— 



2m Q 



(21) 
(22) 



of the equation themselves rather than directly e n . This will be assumed later when 
we study 1/mQ corrections to the static limit. (The explicit factor 1— j3 2 /mq can be 
eliminated by properly rescaling t.) 

The static limit is obtained neglecting all terms suppressed by inverse powers of 
™q 1, 



T / x m 2 -(3 2 / \ t T . N /3 2 
e n Mt) = — 2T"*n(*) + 2^ n{t) - 2 



00 , tyJs) 
as 



(t- 8 y 



(23) 



This is a stationary Schrodinger-type equation for the one-dimensional system with 
the static Hamiltonian 



7~l = 7Y S 



m 2 -P 2 t 



2t 



+ 2 " 



P 2 



ds 



2 Jo (t-s) 2 



(24) 



Similar to Eq. (|6]), the expectation value of the Hamiltonian over a state \l/ n can be 
written in the form 



(n|Wstat|n) 



dt 



m 2 t 
~2t + 2 



I? 



ds dt 



(tt w (s) - %JM 
(s-t) 2 



(25) 



2 Extcnding the interval to [0, oo) introduces errors that are only 0(l/mg). 



5 



Since + | > m, Eq. ( p5|) suggests that e n > m. This lower bound follows also 
directly from the general lower bound Eq. (|7|). 

The large-t asymptotics of the static eigenfunctions is obtained directly: 



p2 F (n) 



where 



_p( n ) 



df&Jt) = lim 

m.Q— +00 



/ ™Q 



(26) 



(27) 



and / n is the usual annihilation constant of a meson ||. The combination FW has 
finite large-iV c and large-mQ limits. Together with various moments (integrals of 
with powers of t) it plays an important role in the heavy quark expansion. For 
future use, it is convenient to define a set of operators acting on ^(t), 



POO 

= / dtt^(t) , Jo = J, 

Jo 



and matrix elements 



F! 



(n) 



dtt l V n (t) . 



(28) 
(29) 



Then, for example, (F^ n A = (n\J \n). Of course, the integrals in Eqs. ( |28| ) and 



literally converge only for — 1 < I < 2. 

The static analogue of the operator K in Eq. (H) takes the form 



00 

ds 



and its commutator with the Hamiltonian 



[H,K] 



m 2 r°° ds 



2t Jo 



1 

2 Jo 



s-t ' 



^ 7712 T 1 T 
as = Ji Jn . 

2t 2 ° 



(30) 



(31) 



For any energy eigenstate n, the equation (n\[H, K]\n) = holds, and we thus find 

2 



or 



m 



m 



dt 

T 

00 dt 



*n(*) 



/ —*„(*) = -P n / dt* n (t) 

Jo t Jo 



(32) 



(33) 



i.e. mF^i = —P n F^ n \ Hereafter we call states for which P n = — 1 and P n — 1 P- 
odd and P-even, respectively. Of course, Eq. (|3~3D can also be obtained directly by 
applying the heavy quark limit to the finite mass parity relation Eq. (|T2|), which is 
the reason why corrections to the static limit do not spoil the parity classification. 
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Another relation among the moments of the wavefunction can be obtained by 
integrating the static equation (E3T) from to oo, yielding 



2e n F^ = m 2 F^ + F[ n) , 
which with the help of the above parity relation can be written as 

Fl n) = (2e n + mP n )F^ . (34) 
Non-diagonal matrix elements of K are also expressed in terms of e and F: 

PnPl-l 1 



{n\K\l) 



p{ n )p( l ) 



2 e n -ei 

Certain useful relations emerge also from the commutation relation 



(35) 



(36) 



and the obvious commutator [t, K] = — J . Thus, for example, [t, [t, H]} = —^~Jo and 

hence (n\[t, [t, H\\\n) = -f (i* n ))M 

In the following, we will derive a tower of relations among the moments of \E^. 
For this purpose, we consider the operators 



A direct computation yields 

n — 1 



t n - 



dt 



dt 



[D n ,H] 



nt n - x U + (m 2 -(3 2 )t n - 2 + 



n+1 



D . 



2 n-3 



(37) 



-t n -yE (k+l)(n-k-2)t M J k . 

(38) 



For n = 0, 1, 2 the last sum is absent. The first few relations take the form 
r, H = 



dt' 



t 2 - K 
t 3 ^ H 

+4 d nj 



m 2 -P 2 , 1 
2? h 2 

-n + t 

-2tH + \(m 2 



3 + 2 



■p 2 ) + it 



n 
n 
n 
n 



-3t 2 U + (m 2 -f3 2 )t + 2t 3 - f J 
-At 3 7i + \{m 2 -(3 2 )t 2 + U A - (3 2 tJ - p 2 Jx n 




1 

2 
3 
4 



(39) 



Taking the expectation values of the operator relations in Eqs. (^) we obtain the 
moments (t n ) in terms of the bound-state energies and decay constants F: 



n 
n 
n 
n 



1 

2 
3 
4 



(t) 
3{t 2 ) 
A(t 3 ) 
5(t 4 ) 



Ait) 2 - (m 2 -{3 2 ) 

6(t 2 )(t)-2(m 2 -/3 2 )(t)+(3 2 F 2 

8{t 3 )(t)-3(m 2 -(3 2 )(t 2 )+A(3 2 F F 1 



(40) 



3 This relation is the so-called fourth sum rule for the Darwin operator. 
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Note that n > 5 does not literally lead to meaningful relations since they would 
involve divergent terms. The case n = yields the relation 



{m-F) J Q f* 2 (t) = 1, (41) 

i.e. the virial theorem Eq. ([19]) for the light quark in the limit where the other quark 
is static. As we have mentioned above, it is literally valid at m > (3, for smaller 
m it can be understood, for example, as an analytic continuation in the mass of 
light antiquark or in a subtracted form. As we will show later, these moments are 
important in the 1/rriQ expansion since the leading I/tuq corrections in Eq. flSip are 
simple powers of t. 

The relation (t) = e n is also a direct consequence of the virial equation (|19"D (the 
one which involves rriq) expanded to the first nontrivial order in 1/itlq. It was first 
derived in Ref. f| (see also [0]). In QCD the bound-state energies e n are usually 
denoted by A n . 

The next important parameter of the heavy quark expansion is the kinetic expec- 
tation value fi^ = 2A ^ h (HQ\Q{iD) 2 Q\Ho) . In the limit uiq — > oo the operator of 

the spacelike momentum takes the simple form, since the Z-graph contributions can 
be neglected, the light-cone combination of momentum is given by xMh q and the 
time component can be excluded using the equation of motion DqQ = rriQ Q. In the 
't Hooft model, in the limit mq — > oo this immediately leads to 

/4 = / df(f-f) 2 * a (f) = (t 2 )-(t) 2 . (42) 
Jo 

Combining the n = 2 relation in Eq. ( |I0|) with the virial equation (n — l) one finds 

3/i 2 = A 2 - m 2 + /3 2 , A = e n , (43) 

for any bound state. 

In the absence of actual chromomagnetic field in D = 2 the next operator is 
represented by the Darwin term 

P 3 d = ^(H Q \Q(-hDE)Q\H Q ) = ^f- (44) 

(the last relation is obtained using the equation of motion for the gauge field, and 
factorization valid at N c — >oo). At the same order a nonlocal zero- momentum cor- 
relator of the kinetic operators appears as well, which will be addressed in the next 
section. 

Before proceeding to the IW functions, let us mention an upper bound on the 
decay constants F^ n \ It is obtained using one of the Sobolev's inequalities bounding 
the Lqo norm via Li and L\ norms in one dimension: 



dz \f(z)\ 2 



1/4 



dz \f'(z)\ 2 



1/4 



for any a. (45) 
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In terms of the Fourier transform of f(z), ^(t) it takes the form 



dtV(t) 



< 7T ' 



dt^ 2 (t) 



iV* r 



| dtt 2 ^ 2 (t) 



1/4 



Applied to the static wavefunction it reads 



" 2 ,2 , ,J2\ " 1 



7T 



4A -m 2 +/^ 



(46) 



or 



1/2 



(47) 



These bounds are aD = 2 counterpart of the bounds discussed in QCD in Ref. ||. 

Since the inequality in Eq. (|46|) is saturated only by functions of the form l/(c-N 2 ), 
a solution of the 't Hooft equation cannot saturate the bound (f|7|). It is interesting, 
however, that for the ground states with light spectator quarks, m ^ (3 the decay 
constant F numerically almost saturates the bound, within only a few percent. 

The operator t plays a special role for the static equation (p3|): the first and the 
last terms in H s t a t are homogeneous functionals of rank —1 with respect to t, while 
the term oc t/2 has rank +1. It breaks the dilatation invariance of the eigenstate 
problem. This operator is the analogue of the operator T> of Ref. || representing the 
part of the full trace of the energy-momentum tensor Q m associated with the light 
degrees of freedom (Sect. II, Eq. (9)). Likewise, there are many relations for various 
observables in the 't Hooft model, involving the operator t. Here we consider matrix 
elements of t. 

Using the commutator Eq. (RO) we write 



{k\t\n) 



1 



e fc — e n 

(3 2 



■(k\[H,t)\n) 



P 2 



4(e fc -e r 



2(e k -e 
(m 2 F^F^ - F (fc) F (n) 



\K\n) 



/? 2 



P 2 



2(e k -e n y 



(k\[H,K]\n) 



\2 - 1 - 1 



;i - p k p n ) , (48) 



where relations (|3TD and (|33| ) have also been used. This can be cast into the form 

(3 2 „(h\ „|V| ( 1 — ( — l) k ~ n 



(k\t-H\n) 



2(efc — e n 



^0 ^0 



(49) 



which embeds both k = n and k^n (we have used the fact that PkP n = (—l) k ~ n ). 
The above equation shows that the operator t— Tt is P-odd, i.e. its matrix elements 
do not vanish only between the states of opposite parity. This makes sense since, in 
the static limit, t — H is simply QiD z Q. 

Using the second of the commutation relations fl39"|) we have for the non-diagonal 
matrix elements of the dilatation operator D 



. d I X 



e 1 



2(e k -e r 



o 



o 



wrong parity 

1 right parity 



(50) 
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These matrix elements determine the so-called oscillator strengths - the Small Veloc- 
ity (SV) transition amplitudes between the heavy quark states (usually the ground 
and the "P-wave" states in actual QCD). Eq. (fffi) allows one to prove an important 
symmetry relation for the IW function. 



3.1 IW function 

The IW function determines the transition amplitudes between two heavy-quark 
states induced by a current bilinear in two heavy quark fields. In the present context 
it can be defined as a diagonal scalar current in the heavy quark limit rriQ — ► oo: 

-^(H^{p')\QQ\H^{p)) = U (WL\ t (gl = (W) s w (51) 
2\/PoPo \ M Hq J M H Q 



In the 't Hooft model the IW functions are given by the following expression in terms 
of the static wavef unctions: 

u = - — 2 — r dt v k (t) ^ n ([w =f Vw 2 - i]t) = 

l + w± vw —1 Jo v ' 

r dt*J-L)*n(^t) , (52) 



1 + z J o \ \/z 

where 

1+z 2 



tr - -- — , z = w ± Vw 2 — 1 • (53) 



The expression for the IW function was obtained in Refs. and |10 



Let us note that each value of w ^ 1 can be represented by two different values 
of z corresponding to two possible values of the square root in Eq. (|52|) , such that 
Ziz 2 = 1. They must yield the same value of £, up to a sign: 

U(z) = P n PkU(l/z) , (54) 

which, for n^k looks like a miraculous property of the 't Hooft equation ||. Alter- 
natively, the above property can be written as 

£nk(w) = P n P k £ kn (w) . (55) 

Now we can demonstrate it explicitly.0 

Using the fact that D = is the generator of scale transformations, i.e. 

e ln(a)D /(i) = e ln(a) ^/(t) = f(at) (56) 



i We thank R. Lebed for its cross checks in the numerical computations. 
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for an arbitrary function /(£), the IW function can be written in the form 



ink{z) 



2y/j 
l + Z 



(k\e 



(£>+§) In* 



D = t- 



dt 



(57) 



Note that the operator -D+| is antihermitean (antisymmetric), so that ||e ( D+ ^ lnz \n) \\ 
= || |n)|| . This ensures the so-called first sum rule expressing the unit probability of the 
transition to arbitrary final state in the heavy quark limit (in the SV approximation 
it is known as the Bjorken sum rule). 

To calculate the diagonal matrix elements of D one can use the identity 



/ dt *'(*) *(t) /(f) = ~ Jdt j t (* 2 (f) f(tj) - * 2 (f) f(t) 



"2 </'(*)) 



valid for arbitrary fit). In particular, it shows that the expectation values of D + | 
vanish. Together with relation (|50|) we see that only even powers of D + | survive in 
the exponent in Eq. ( 57f) when \n) and |/c) have the same parity, and only odd powers 
contribute if the parity of the two states is opposite: 
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nk 



2y/i f (jfe|cosh[(D + i)ln«]|n) 
1+2 1 (fe|Binh[(D + |)lnz]|7i) 



Ti- 
ll - 



■k 
■k 



even 
odd 



(59) 



This proves the symmetry properties Eqs. (|54|) , ( 55|) and ensures that the IW func- 
tions are analytic at vv' = 1, in spite of the branch point in z as a function of vv' . 



3.2 SV sum rules 

Important constraints on the transition amplitudes between heavy flavor hadrons 
and on the parameters of the heavy quark expansion follow from the sum rules, in 
particular, in the small velocity (SV) limit. In 1+1 dimensions the first four sum 
rules in the heavy quark limit take the form 

2 _ ' V"^ 2 

Pk ~ T _ / j T lk 
4 I 

\^ = Ete- 6 *)^ ( 6 °) 



k 



p%) k = E(^-^)V fc - (si) 
i 



The so-called 'oscillator strengths' r parameterize the transition amplitudes into the 
opposite-parity states in the SV limit, and p 2 denotes the slope of the elastic IW 
function: 

2^W7,.0I*> = n k e,„v», —( k (v)\Q lo Q\k(0)) = l-fgL + ofr) (62) 
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(v is the velocity of the final state hadron). Therefore, 

Tlk = {llt ft + l lkh pl ~l = {kl { t ft + l) lk) - (63) 

The first sum rule then becomes obvious being a consequence of completeness of the 
eigenstates. Other sum rules are straightforward as well. 

Consider, for example the third sum rule for the kinetic operator. Using the 
commutator with n = 1 in Eq. (^), we have 

E4(Q-e,) 2 = J2(k\t-n\i)(i\t-n\k) = (k\e\k) - mm? = (fi) . (64) 
i i 

Similarly, we get for the second, "optical" sum rule 

= -12(k\t\l)(l\tf t + l\k) = = f • (65) 

The fourth sum rule for the Darwin operator can be directly obtained by inserting 
the complete set of states into the commutator [t, [t,7~t]] = —\Jq '■ 

E*»te- c *) 8 = ~(k\[t,[t,n\\\k) = ^ F p 2 = ( P i) k . (66) 



4 l/ m <3 expansion 

In practice, it is often necessary to account for the first few 1/ttiq corrections to the 
static limit tuq oo, since in actuality these effects are non- negligible not only for 
charm, but even for beauty hadrons. In studies of the 't Hooft model there appears 
an additional motivation: the available numerical approaches often apply only to the 
finite quark masses. The solution of the static equation (|23|) is approximated by the 
solution of the generic finite-mass 't Hooft equation where mg is taken large but 
finite. For computational reasons mg cannot be taken too large, and control over the 
'spurious' 1/rng corrections becomes mandatory even for studies of the pure static 
case. 

In this section we will study the leading I/tuq corrections to the axial decay con- 
stant, meson masses and the kinetic energy of the heavy quark. The I/tuq expansion 
is carried out by applying to Eq. (^Tj) the standard formalism of non-covariant time- 
independent perturbation theory used in QM. Since the leading 1/rriQ corrections in 
Eq. (pT|) involve only powers of t, it is possible to derive exact expressions for the 
first few terms in the 1/rriQ expansion for these observables that depend only on the 
moments of the static wavefunction and of the structure function. Using the results 
from Sect. 3, one can then express the coefficients in the 1/mg expansion solely in 
terms of the static binding energy and the decay constant. 



12 



Although we will perform the 1/ttt.q expansion using old-fashioned time-ordered 
perturbation theory, we will introduce here the corresponding notations which resem- 
ble those used in field theory, where 1/toq corrections to various expectation values 
are given by correlators of the type 

/oo 
dTiT{A(0),5H(r)}\k) . (67) 
-oo 

Heisenberg operators 0(t) are understood as e lHr O(0) e~ lHT ] we assume that the 
Schrodinger operators we deal with, do not depend on r explicitly. 
We then denote for the stationary problem 



/oo 
dr iT{A(0),B(r)}\k) 
-oo 



£ (k\A\n)(n\B\k) + (k\B\n){n\A\k) ^ (gg) 

The similar expectation value can be defined for the time-ordered product of arbitrary 
number of operators. We consider 

H ^H + aA + f3B + -fC + ... 

and put 



( W , B , C ,...}|^-(|;^-.*A7,..) 



(69) 



/ a =/3= 7 =...=0 

Two basic relations hold for such T-products of two operators: 

(k\iT{[H,A],B}\k) = (k\[A,B]\k) (70) 

and 

(k\iT{AH, B}\k) = e k (k\iT{A, B}\k) - (k\AB\k) + (k\A\k)(k\B\k) , 
(k\iT{HA,B}\k) = e k (k\iT{A,B}\k) - (k\BA\k) + (k\A\k) (k\B\k) . (71) 

The first relation has a transparent meaning: since [H, A] = —i^rr, Eq. ([7(]) is a form 
of integrating by parts: 

(*l /drT |^,S(0)J|*> = (fc| /dr A T{A(r))5( o) } |fc) + (A;|[A(0),i?(0)]|A ; ). 

(72) 

The obvious relation 

(A;|zT{^,A}|A;) = 
holds as well, which will be used later. 
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In the static limit, (t) equals to the bound-state energy e. It is often necessary 
to know how (t) changes under various perturbations 57i. The answer is readily 
obtained using Eqs. fl39|) and (f70|): 



6t kk = (k\iT{t,SH}\k) = -(k\[t^,5H]\k) . (73) 
Since D\ = t4i generates scale transformations in t, one finds 

[t^O] = Dim[0]-0. (74) 

Say, [t-^,t l ] —It 1 . Therefore, any perturbation which is a homogeneous rank-/ func- 
tional of t, satisfies Jl^ 

5(t) = -l(6Hi) ■ (75) 

For example, for 57i = Xt one finds 5(t) = —A, a result that can be easily verified by 
direct evaluation, since the exact result reads (t)\ = (t)/yl + 2X. 

The same property holds for non-local operators as well. For example, 

(k\iT{t,A,B}\k) = ~{D[A]+D[B] + 1) (k\iT{A, B}\k) , (76) 

or 

(k\iT{(t-2H),A,B}\k) = -(D[A] + D[B]-1) (k\iT{A, B}\k) , 

where D[B] denote the t-dimension of operators A and B. 

The above properties parallel the relations for the operator of the trace of the 
energy-momentum tensor for the light degrees of freedom in actual QCD discussed 
in Ref. |J. This similarity was elucidated in the previous section. 

The mQ-suppressed terms in the r.h.s. of Eq. ( pID playing the role of the pertur- 
bation 57i are given by t 2 /2rriQ, t 3 /2mg, etc. Therefore, in the 1/toq expansion 
one typically needs to compute T-products Eq. ( j68|) with the operator t 2 (or t 3 , in 
higher order). As exemplified above, this can be done directly using the relations in 
Eqs.QTOD and Qnj): 

(k\zT{t 2 , A}\k) = ~(k\[D 2 , A]\k)-~e h (k\[D u A]\k)-~(k\tA\k)+^e k (k\A\k) , (77) 

(k\iT{e-^,A}\k) = ~(k\[D 3 , A]|A;)- efe (fc|[ J D 2 ,A]|A;)-i(4e2-m 2 +/3 2 )(A;|[ J D 1 ,A]|fc) 
- ^(k\t 2 A\k)-2e k (k\tA\k) + l -(%e 2 k -m 2 +(3 2 )(k\A\k) , (78) 

etc. 

As an application of these relations, we obtain 
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(k\iT{t 2 ,t 2 }\k) = --(2(t 3 ) + (t)(t 2 ) 



6 44-fe k (m 2 -p 2 ) + 2 -p 2 F 2 



(79) 
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This correlator governs the 1/tuq corrections to the average (t 2 ) which, in turn, 
determines the kinetic expectation value in the static limit. Likewise 

(k\iT{t\t l }\k) = - 2 l±±(k\t l+1 \k) - ^(l-l)e k (k\t l \k) . (80) 

Similarly we get the analytic expression for the 1/mq corrections to the decay 
constants F^. Indeed, (F^) 2 = (k\J \k), and 

si /mQ (F {k} ) 2 = ^(m{t 2 ,j Q }\k) = 

1 -(k\ - 2 -[D 2l J ]\k) + ~ {t k (k\iT{t, J }\k) - (k\tJ \k) + (k\t\k)(k\J \k)} = 



2rriQ 3 2niQ 3 

_1_ {-^F^-I^W,, JU|*>-f fftf'+l^ff)'} = -_j_|tf>*<» 

(81) 

where we have used that [Di, J ] = J and [D 2 , Jo] — 2Ji- Thus, 



5 1/mQ F^ _ 1 2^ 2(2e k +mP k ) 



F( fc ) m Q 3 3m c 



(82) 



In the last equation we used relation Eq. (|34D for Fi"\ The 1/rriQ corrections to F 
turn out significant (c.f. Refs. [|, |Tl|1 ). 

It is often advantageous to define the axial decay constant via the pseudoscalar 
current rather than the axial current: 

mi 5 Q\k) = (83) 



2M Hq " ^ 2 
then 

f ~(k) = fW _M ]L _ ^ p{k) = F(k) _J4_ 

rriQ+m niQ+m 

Similar to what is observed in actual QCD, the I/tuq corrections to the ground-state 
F are smaller, 

fymg F {k) e k +m{3 + 2P k ) 

W) ~ 3m Q • m 

The analytic expression ( |8"2"D agrees with the numerical computations performed 
in Ref. H for the ground state. In terms of c k introduced there to quantify these 
preasymptotic corrections, 

5 

Ck = ~£k + mP k . (86) 
o 

In the nonrelativistic case e k — > m holds, and for the negative-parity ground state 
one has the correct limit c = — m/6 |J. For the first excitation, however, one would 
have Ci — > llm/6. 
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It is not difficult to derive the expression for the 1 / wiq correction to the light-cone 



wavefunction itself generated by perturbation 



2mq ' 
1 -(1., 



= *?(*) + — Wif)- (87) 

m Q 



Using the commutators in Eqs. we get 



f <X)U\ - 1 L-2 d ■ » ■ d 



*™ = 3i t di + 2£ ' t di + 2£ 'J*" W ' 

Note that simply replacing wlq by ~ tuq + e when passing from ifk{x) to 

^(t), would amount to adding only ^efc + |J ^(t). The remaining part 

3m^ (t 2 &t ~~ efe *cl + ¥) ac t ua Uy changes the shape. 

Using similar techniques, it is straightforward to obtain explicit I/tuq expansion 
of the meson mass in the 't Hooft model; we consider it through order /3 3 /mg corre- 
sponding to the order discussed in case of QCD ||. A straightforward evaluation of 
the expectation value of Eq. (^) yields (see also fLOf ) 

M _ (t , , (t 2 )oo-(t)l-(3 2 , 4(t 3 ) + 4(t) 3 + (*T{t 2 ,t 2 })-4(t)(t 2 ) ( {3*\ 
M HQ -m Q - (tU+ — + ^ + °U&J" 

(89) 

Here all expectation values correspond to the static limit w,q — ► oo, i.e. do not 
implicitly depend on uiq. Using the previously derived relations, especially Eqs. (|40|) 
and (|79|) , all terms can be expressed only in terms of the static binding energy A and 
the axial decay constant of the respective state. 

Finally, the object usually appearing in the l/mq expansion of the diagonal matrix 
elements to this order, is the zero- momentum correlator of operators 7F 2 = Q(iD ) 2 Q 
which represent the 1 / mg piece of the Hamiltonian, 

-PL = \{k\iT{^}\k). (90) 

In particular, it determines the 1/m.Q variation of the kinetic expectation value 
(k\Q(iD ) 2 Q\k) itself in the actual fmite-mQ hadron. The expression for p 3 ^ in 
the 't Hooft model is most simply obtained using the above mentioned relation 
Q(—iD z )Q = t—Ti which holds for the zero-momentum matrix elements in the static 
limit. Therefore, we simply need to compute (k\iT{(t— H) 2 , (t— H) 2 }\k). In this way 
we obtain 

-2/4 = (k\iT^ 2 ^ 2 }\k) = A{e)+A{tf + {iT{t 2 ,t 2 })-A{t){t 2 )-^-F 2 = (91) 
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-(l<? k --^-f) + \^) 
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This correlator is numerically large. 

We can compare the 1/rriQ expansions discussed above with the general operator 
expansion Ref. R |7| valid in arbitrary gauge theory: 



M HQ -m Q = A 4 



2mn 



--(zT{7f 2 ,7f 2 }) + ^ 

8 4 



B 2 F 2 

+ ..., p 3 D = ^.(92) 



Similarly - and even a bit simpler - is to consider the expansion of the scalar density 
which is given precisely by |7| 



(QQ) 



m Q A 



2M Hq M Hq V 

and, on the other hand, use the similar OPE for (QQ) : 

1 = (Q(iD ) 2 Q)-f3 2 p% 



(93) 



2M, 



\Hq\QQ\Hq) 



1 



2m^ 



my 



(94) 



The explicit computations show that these equations are satisfied with given by 
Eq.fT]). 

It is interesting that it is possible to derive a closed expression for the expectation 
value of the 'kinetic' operator Q(iD z ) 2 Q in terms of the 't Hooft wavefunction for 
arbitrary mass mg: 



2M 



■(H q \Q(iD z ) 2 Q\H q ) 



2M, 



Hq 



x 4 \mQ{mQ — mP) Jo 



dx <p(x) 



(95) 



(this expression assumes a certain ultraviolet regularization of the operator, see be- 
low), where P is parity of Hq. The idea is the following. 

In the rest frame of the meson the expectation value of Q[(iDq) 2 + (iD z ) 2 ]Q is 
simply 2(Q{iDJ) 2 Q) and, therefore 

^—(H Q \Q[(iD ) 2 + (iD z ) 2 }Q\H Q ) =m Q M HQ f dxx V \x) . (96) 

The complementary combination of momentum operators Q[(iDo) 2 — (iD z ) 2 ]Q can 
be determined using the general identity 



d D x QQ{x) 



d D x { g 7o Q + Q 



JiD -m Q ) 2 - (iDf + \o» v G 



2mi 



Q (97) 



valid in arbitrary dimension. In D = 2 one has ^a^G^ = Gi'j 5 with G = \t^ v G^ v ■ 
GVji- Since (QiD z Q) = (QiD iD z Q) = and the light-cone combination (Q(iD - 
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iD z )Q) is again directly expressed via the 't Hooft wavefunction, Eq. ( 97f) yields the 
necessary equation for (Q[(iD ) 2 —(iD z ) 2 ]Q) involving, however the expectation value 
(QGi^Q). 

The gluonic field strength Got can be explicitly written via quark current in the 
light-cone gauge: 

G i = -9lj-J + (98) 

which leads to a non-local four-fermion operator. We note, however, that the current 
J + includes not only the spectator quark q, but also Q"i + ^Q- This term described 
by the bare loop leads to the ultraviolet divergent expression. We simply discard this 
contribution in J + in Eq. (0), which fixes a certain renormalization procedure. Then 
we get 



1 4 Mh q Jo Jo 



1 r l 





(99) 

This finally yields Eq. (|95|) where we have used Eqs. (|14D to express the last term via 
the decay constant of the meson. 

A note of caution must be voiced regarding this derivation, however. In the 
way described above we obtain the bare operator Q{iD z ) 2 Q. It does include a finite 
contribution from the domain of momenta of order ijiq even in the leading order 
in rriQ. On the contrary, in the heavy quark expansion we are interested only in 
the physics originating from momenta essentially below the scale of the heavy quark 
mass itself. The expressions Eqs. (0), (|91| ) refer just to such low-energy effective 
operator. Therefore, in general the literal comparison of the two expectation values 
is not too instructive. It is easy to check that already to the leading order in mg the 
two expressions differ by the amount %r attributed to the domain of momenta ~ niq. 

Here we note an interesting feature of the exact expectation value of the local 
quark-gluon operator QGi^^Q in Eq. (p9|). At mq — > it has an I/uiq singularity 
regardless of the mass of the second quark in the meson: Jq 1 dx/x <p(x) oc 1/wiq at 
mg (3. Yet we know that no appropriate massless physical states exist in the 
model at m q ^ (the ground QQ state has negative parity, and their pairs are 1/N C 
suppressed). The singularity technically emerges due to massless gluon propagator, 
however gluon is absent from the physical spectrum in D = 2.1 

It must be noted, however, that careful treatment of passing to the light-cone co- 
ordinates in the computations of the similar vacuum expectation value (O^Cr^^lO) 
in the 't Hooft model yielded additional terms which canceled the pole observed 



5 We are grateful to A. Vainshtein for informing us of existing examples of similar IR singular- 
ities in physical amplitudes in the absence of contributing massless particles, in low dimensions. 
Reportedly, such a situation is excluded in D > 3. 
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in Ref. Jl2| and led to a finite result at m^ — >• 0. A possibility of similar subtleties 
in the computation of the meson expectation values deserves further studies. We are 
grateful to A. Zhitnitsky for pointing out and discussing this problem. 



4.1 



I/tUq correction at zero recoil 



As another application of the 1/mg expansion, we briefly consider here the second- 
order nonperturbative corrections to the zero-recoil B — > transition amplitude. 
At this kinematic point the deviation from the elastic IW function (which is unity 
here) appears at the level l/m 2 cbl which provides a method of extracting \V c b\. The 
corrections, however, are shaped by strong dynamics at the typical hadronic scale 
and at present cannot be evaluated from the first principles. The existing estimates, 
in particular for the axial B — > D* amplitude, rely on the sum rules derived in 
Refs. @,|TJ: 



D* 



6 



/4 

3m^ 



1 1 

^7 + ^7 + 



mt 



mi 3m c mb 



O 



m f 



(100) 

where are the zero-recoil transition amplitudes to the excited states, oc I/tuq, 
£a is the short-distance renormalization factor, and fi 2 , are expectation values of 
the kinetic and chromomagnetic operators, respectively. One then has [13] 



D* 



A/2 



j4_ 

6m? 



+ 



+ 



+ 



mt 



mt 3m c mi, 



mi 



(101) 



where a positive quantity \ parameterizes the magnitude of the sum of the excitation 
probabilities in the l.h.s. of the sum rules, in terms of the local operator term in the 
r.h.s. of Eq. (]T00|) : 



E 1^1 

k^O 



x- 



/4 

3m? 



+ 



mi ~* "' 2 



mi 3m c m}j 



(102) 



The expressions for the excitation amplitudes were elaborated in Ref. ||. Following 
Ref. [13 1, existing numerical estimates of Fd* assume (somewhat arbitrarily) that x 
can vary up to 1, that is, x — 0-5 ± 0.5. 

We computed x analytically in the 't Hooft model. Since spin and chromomagnetic 
field are absent in two dimensions, only the kinetic operator acts here. The sum rule 
takes the form 



F*+Y\F k \ 2 



k^0 



u - 



1 

m a 




m h j 



(103) 



and the excitation amplitudes Ff, to the leading order are given by 

1/1 1 \ (k\Q^ z Q\0) 



m r 



m b 



(104) 
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Similarly, x is defined through 



El^l 2 = X-(^-^) 2 ^. (105) 



The sum of F% can be computed using the same technique as was elaborated in 
the previous sections. However, the corrections to the amplitudes we consider are not 
expressed anymore via only positive integer moments, and include expectation values 
of operators with higher derivatives. Yet they can be expressed in terms of the slope 
of the IW function p 2 . We give here the final result 

10 5 e 2 n 4 / „ 1\ 5 5 m 2 -(3 2 4 / „ 1\ , , 

For light spectator quark the value of x turns out to be about 0.55. This is 
surprisingly close to the central value guestimated in Ref. [O] for the case of actual 
QCD. 



5 Summary and Outlook 

We have studied the 't Hooft model in the mq — > 00 limit. Our main result are 
exact relations for the heavy quark kinetic energy, as well as moments of the 't Hooft 
wavefunction and of the structure function in this limit, which allow to express these 
observables in terms of only the heavy quark binding energy A and the axial decay 
constant F k . 

In the 't Hooft model, these moments appear as coefficients in the 1/toq expan- 
sion for various observables. As an example, we calculated coefficients in the I/tuq 
expansion of meson masses, decay constants and heavy quark kinetic energies. Using 
the above relations, we were able to express the corresponding 1/tuq coefficients in 
terms of A and F k . 

Likewise, we derived the expressions for the oscillator strengths and verified a set 
of the SV sum rules in the heavy quark limit. As an application of the developed 
1/toq expansion, we computed the nonperturbative 1/rriQ corrections to the zero 
recoil B — » D transition amplitude. 

Although the 't Hooft model is in principle "numerically solvable" , many observ- 
ables can only be determined with very limited precision in practical calculations. 
This is particularly the case for observables in the limit where one of the quarks 
becomes heavy. In this case the 't Hooft wavefunction becomes extremely asymmet- 
ric and many numerical techniques, which are otherwise rather powerful for finite 
quark masses, fail to produce numerically reliable results. For this regime, where the 
heavy quarks are not infinitely heavy, it is often advantageous to perform the I/tuq 
expansion beyond the leading order. The corresponding expansion coefficients that 
we derived involve only properties of the 't Hooft wavefunctions in the static limit. 
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Moreover, through the use of exact relations, the expansion coefficients can be ex- 
pressed in terms of mq — > oo properties of the wavefunctions that can be calculated 
numerically with sufficiently high accuracy. 

The developed analytic 1/mg expansion allows to carry out precision studies, in 
the framework of the 't Hooft model, of such an intriguing and poorly understood 
phenomenon as violation of local duality in heavy flavor decays |14|]. The question 



of its magnitude has a particular phenomenological significance in the domain of 
moderately heavy quarks, where reliability of the asymptotic expansions is unknown 
a priori, and numerical computations are unavoidable. 

The developed technique can be used to test, on the example of the 't Hooft 
model, various approximations routinely applied to the actual beauty decays. Some 
of them will be reported in Ref. [^4[ , while others further deserve dedicated studies. 
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